Polyakov loops and spectral properties of the staggered Dirac operator 
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We study the spectrum of the staggered Dirac operator in SU(2) gauge fields close to the free 
limit, for both the fundamental and the adjoint representation. Numerically we find a character- 
istic cluster structure with spacings of adjacent levels separating into three scales. We derive an 
analytical formula which explains the emergence of these different spectral scales. The behavior on 
the two coarser scales is determined by the lattice geometry and the Polyakov loops, respectively. 
Furthermore, we analyze the spectral statistics on all three scales, comparing to predictions from 
random matrix theory. 
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I. INTRODUCTION 

One of the goals of quantum chromodynamics (QCD), 
the fundamental theory of quarks and gluons, is the cal- 
culation of the hadronic mass spectrum from first princi- 
ples. For this purpose, a nonperturbative regularization 
of QCD can be formulated on a space-time lattice, which 
makes the theory mathematically well defined and ac- 
cessible for numerical simulations. In order to describe 
fermionic fields, one needs a discretizcd Dirac operator, 
which can be constructed in different ways. In particu- 
lar, the staggered or Kogut-Susskind Dirac operator [l| is 
widely used as it is computationally cheaper than other 
options. (We shall not enter the debate of the rooting 
issue [1, 3 here.) 

During the past 15 years it has been shown that some 
universal properties of the QCD Dirac spectrum can be 
described by a version of random matrix theory (RMT) 
[J,i5[] which incorporates the chiral symmetry of the mass- 
less Dirac operator and is accordingly called chiral ran- 
dom matrix theory (chRMT) [I, 0]. In chRMT, one 
models the Dirac operator by a block off-diagonal matrix 
with random entries, respecting the global symmetries of 
the massless Dirac operator. The path integral over the 
gauge fields is then replaced by averaging over an ensem- 
ble of such matrices. Chiral RMT reproduces low-energy 
sum rules Q and yields accurate predictions for, e.g., 
the microscopic spectral density Q, the distribution of 
the smallest eigenvalues [10, [U , level spacing distribu- 
tions, and other short-range spectral correlations. Pre- 
dictions from chRMT have been successfully compared 
with numerical results from latti ce gaug e theory in many 
different settings [H, M, 113, M, [M [II E [H, S, [Si • 

When using RMT to predict spectral properties of 
complex quantum systems, it is essential that the ran- 
dom matrices have the same antiunitary syrnmctries as 
the system to be modeled, see, e.g., Ref. [22]. One dis- 
tinguishes different symmetry classes and corresponding 
ensembles of random matrices, labeled by their Dyson 
index Pd = 1, 2, or 4 d]. 

The staggered Dirac operator on the lattice exhibits 



the peculiar feature that its symmetry properties can be 
different from that of the continuum theory. In partic- 
ular, for gauge group SU(2) with fermions in the funda- 
mental representation and for gauge group SU(A^) with 
fermions in the adjoint representation, the (3d = 1 and 
/?_D = 4 cases are interchanged compared to the contin- 
uum Dirac operator (see below). This implies that the 
spectral properties of the staggered operator are differ- 
ent from those of the continuum Dirac operator. Ac- 
cordingly, a transition is expected to take place in the 
continuum limit. The first indication of such a transition 
has been reported in Ref. [23J . 

In this work, we study a related but different case, 
namely, the staggered Dirac operator for fermions in the 
fundamental and adjoint representation of SU(2) in the 
free limit. This limit is approached by increasing the 
Wilson gauge action parameter (3 = A/g^ at fixed (or 
mildly varying) lattice size, i.e., the lattice spacing, and 
thus the physical volume, shrinks to zero. 

In the free limit the dynamics becomes integrable and 
therefore one no longer expects RMT statistics. Conse- 
quently, in this limit we do not expect the same transi- 
tion between RMT symmetry classes as in the continuum 
limit described above. Instead, generically, the eigenval- 
ues should be uncorrelated, like numbers drawn from a 
Poisson process ;^24i, ,25i]- It turns out that in the case 
we are studying, the situation is somewhat more compli- 
cated: Close to the free limit it is possible to disentangle 
three different scales that appear in the separations be- 
tween the eigenvalues. 

First, the eigenvalues form well-separated plateaux 
centered at the eigenvalues of the free staggered oper- 
ator. 

Second, we observe an internal structure of the pla- 
teaux: The eigenvalues arrange themselves in clusters of 
eight eigenvalues each. We will show that, for each con- 
figuration, the location of these clusters can be predicted 
from the knowledge of only four real variables, i.e., the 
averaged traced Polyakov loops in the four lattice direc- 
tions. 

Third, on the finest scale, the eigenvalue fluctuations 
within clusters can be described in terms of chRMT. We 



will present numerical results for the level spacing densi- 
ties which agree with the chRMT prediction for all values 
of the Wilson gauge action parameter /?. 

This article is organized as follows. In Sec.[n]we briefly 
introduce some basic notions of random matrix theory 
and discuss the different antiunitary symmetries which 
are relevant for the Dirac operator in the continuum and 
for the staggered Dirac operator on the lattice. In Sec. lIIII 
we derive an analytical prediction for the staggered Dirac 
eigenvalues in certain gauge field configurations close to 
the free limit. Section IIVI begins with a comparison of 
numerical data from lattice simulations to our analyti- 
cal prediction and continues with an analysis of spectral 
statistics on all three scales; i.e., we study the distribu- 
tion of level spacings within clusters, between clusters, 
and between plateaux. Wherever appropriate we com- 
pare to the predictions from chRMT and from a Poisson 
process. Throughout Secs.HIlto lIVI we discuss the funda- 
mental and adjoint representation side by side. In Sec. V 
we conclude with a discussion of the implications of our 
findings. Preliminary results of this study have been pre- 
sented in Ref. HI. 



However, the Dirac operator enjoys invariance with re- 
spect to an antiunitary transformation if the fermions 
are either in the fundamental representation of the gauge 
group SU(2), or in the adjoint representation of SU(A^) 
with N > 2 arbitrary. These two cases are discussed in 
more detail in the following. 



A. Fundamental representation 

In a nutshell, the antiunitary invariance of the Dirac 
operator with SU(2) gauge fields and fermions in the fun- 
damental representation, i.e., QCD with two colors, is 
based on the fact that the generators are Ta/2, and the 
Pauli matrices r^ possess the following complex conjuga- 
tion property: 



-T2TaT2 



(2.3) 



The antiunitary symmetry operator, however, is realized 
in different ways in the continuum and in the lattice for- 
mulation with staggered fermions. 



II. ANTIUNITARY SYMMETRIES AND RMT 
ENSEMBLES 

The notion of universality, as commonly used in the 
context of RMT and the analysis of spectra of complex 
quantum systems, means that spectral statistics can be 
described by an appropriate ensemble of random matri- 
ces, which shares the symmetries of the system under 
consideration. Depending on the presence of antiunitary 
symmetries, the entries of matrices of the ensemble have 
to be either real, complex, or quaternion real. The asso- 
ciated Dyson indices are (3d = 1, 2, or 4, respectively. 

In particular, RMT yields a prediction for the univer- 
sal quantity P{s), the probability density for the unfolded 
nearest-neighbor spacings s (see Scc. lIVCTj for a discus- 
sion of unfolding) . This prediction is well approximated 
by the Wigner surmise. 



P{s) 



asP-e-'^ 



where 



a = 2 



r^°+i (/3d/2 + 1) 



b = 



r^ {I3nl2 + 1) 



rfo+2 ((/3^ + i)/2) ' r2 {{(3d + l)/2) 



(2.1) 



(2.2) 



The Wigner surmise is the level spacing density for en- 
sembles of 2 X 2 matrices with Dyson index (3d- 

For the QCD Dirac operator D the RMT description is 
formulated in terms of matrices which reflect the chiral, 
flavor, and antiunitary symmetries of D [27|. For gauge 
group SU(A^) with iV > 3 colors and fermions in the 
fundamental representation, D generically has complex 
elements and does not commute with any antiunitary 
operator. Accordingly, universal spectral correlations are 
described by the chiral unitary ensemble (chUE), labeled 
by the Dyson index (Jd = 2. 



1. Continuum 

In the continuum, the anti-Hermitian massless Dirac 
operator is defined as 



D = 7m^m = 7m {^^. + ^A^Ta) 



(2.4) 



with Ta = Ta/2 for the fundamental representation of 
SU(2). It anticommutes with the chirality operator 75, 
and therefore its nonzero eigenvalues, which are purely 
imaginary, come in complex conjugate pairs. 

Using Eq. (|2.3|1 one easily verifies that the operator D 
is invariant under an antiunitary symmetry. 



[C75T2i^,Z?]=0, 



(2.5) 



where C = 7274 is the charge conjugation matrix, and K 
denotes complex conjugation (in the position representa- 
tion). Note that C75 acts on the spinor indices, whereas 
T2 acts in color space. Since {C"f^T2KY ~ 1, it follows 
that D can be made real by a basis transformation that 
does not depend on the gauge configuration [23, 1281 . Ac- 
cordingly, the RMT description for two-color QCD and 
fundamental fermions is formulated in terms of the chiral 
orthogonal ensemble (chOE), characterized by (3d ~ 1- 

2. Lattice 

The staggered or Kogut-Susskind Dirac operator for a 
hypercubic lattice of finite spacing a in d dimensions is 
given by 



1 

(2.6) 



M=l 



with r?^(x) = (_l)i:,.<M^- and Uf,{x) G SU(2). On the 

lattice the operator S = Sx,yi—i)^'^=^'^" plays the same 
role as 75 does in the continuum: Since {DkSj'S'} = 0, 
the eigenvalues of Dks also come in (purely imaginary) 
complex conjugate pairs. 

The staggered Dirac operator, which is widely used in 
numerical simulations (because it maintains a remnant of 
chiral symmetry, partially solves the doubling problem, 
and is computationally cheaper than other lattice Dirac 
operators) exhibits the peculiar feature that its antiuni- 
tary symmetries are different from those of the continuum 
Dirac operator [H, [H, [s^l ■ Because of Eq. dM]) the links 
obey Ufx{x) = T2U*{x)t2. Since the 7 matrices have been 
replaced by real numbers, rj^(x), no charge conjugation 
is required in order to compensate for the complex conju- 
gation. Therefore, I?ks is invariant under the following 
antiunitary symmetry: 



[T2K, D^s] = 



(2.7) 



As (t2A')^ ~ —1, it follows that Dks can always be writ- 
ten as a quaternion real matrix [22| . Hence, the chiral 
symplectic ensemble (chSE), with fio = 4, is used to de- 
scribe its universal properties. Another consequence of 
invariance with respect to an antiunitary transformation 
with square —1 is Kramers' degeneracy, see, e.g., 22, 31[; 
i.e., all eigenvalues have (at least) multiplicity two. This 
degeneracy has to be removed by hand before one dis- 
cusses spectral correlations. 



B. Adjoint representation 

In this case, the antiunitary symmetries are deter- 
mined by the purely imaginary nature of the matrix ele- 
ments of the generators in the adjoint representation (be- 
ing the structure constants) of the gauge group SU(iV). 
As in the previous case, the way this symmetry is realized 
in the continuum and for the staggered Dirac operator on 
the lattice is different. 



1. Continuum 

The purely imaginary adjoint generators induce a real 
covariant derivative in the continuum Dirac operator of 
Eq. (|2.4p . which is therefore invariant under the following 
antiunitary symmetry: 



[C75i^,D] = 



(2.8) 



Since {C^c^KY = ~1j it is possible to recast D into real 
quaternionic form. This implies that the Dirac spectrum 
can be described in terms of the chiral symplectic ensem- 
ble (chSE), labeled by 13d = 4. 



2. Lattice 

For fcrmions in the adjoint representation, the stag- 
gered Dirac operator Z?ks is explicitly real, because such 
are the Ufj,{x) link matrices appearing on the right-hand 
side of Eq. (|2.6p . which now take values in the adjoint 
representation of SU(2). This implies that the appro- 
priate RMT ensemble is the cliOE, characterized by the 
Dyson index /3d = 1. 



III. SPECTRAL PREDICTIONS FROM 
POLYAKOV LOOPS 

The spectral properties of the Dirac operator are rel- 
evant to the two major unsolved nonpcrturbative prob- 
lems in QCD: confinement and chiral symmetry break- 
ing. On the one hand, the eigenvalue density is related to 
the chiral condensate, the order parameter for the chiral 
phase transition, by the Banks-Cashcr relation [32|. On 
the other hand, the average value of the Polyakov loop, 
which, in the quenched case, is an order parameter for the 
confinement-deconfinement transition, can be expressed 
through sums of the eigenvalues of the lattice Dirac oper- 
ator with different boundary conditions [33|, |3J, |35|, |36| . 

On commensurate lattices, i.e., on lattices for which 
the numbers ip/2, ^ = l,...,d, are rationally depen- 
dent (the staggered Dirac operator only makes sense on 
lattices with even L^), the spectrum of I?ks in the triv- 
ial vacuum, i.e., for the field configuration with all [/^(.t) 
equal to unity, is highly degenerate. Our analysis of 
the spectrum close to the free limit shows that the way 
in which this degeneracy is (partially) lifted can be ex- 
pressed in terms of the traced and averaged Polyakov 
loops in all directions. 



Pu = - tr 
^ 2 



Y[ U(x + njl) 



(3.1) 



where the x average is over the whole lattice. For each 
configuration, this effect can be predicted by calculating 
the spectrum obtained from a vacuum configuration with 
uniform link variables in each direction, taking values in 
an Abelian subgroup of the gauge group, and yielding 
the same averaged Polyakov loops as the original config- 
uration. 



A. Fundamental representation 

For the trivial gauge vacuum (all links set to 1 or any 
jauge transform thereof), the eigenvalues of the stag- 
gered Dirac operator read 



±i 



\ 



d 

E 



sm 



27r 



(fcp 



(3.2) 



with the wave numbers taking integer values < fc^ < 
L^/2 and c^ = (c^ = 5) for (anti)periodic boundary 
conditions for the fermionic wave function in direction fi. 

However, the system also admits vacua in different cen- 
ter sectors, which can be labeled by the traced Polyakov 
loops Pp = ±1 [for SU(2)]. The sign of the latter can 
always be inverted by a multiplication of all the links in 
the fi direction in a given fixed-x^ slab by —1. In the 
Dirac operator this can be compensated for by switch- 
ing from periodic to antiperiodic boundary conditions 
(or vice versa) in the direction fi. Therefore, when we 
consider a vacuum where the Polyakov loop in direction 
fi is —1, we can equivalently set c^ = ^ (c^ = 0) for 
(anti)periodic boundary conditions of the Dirac operator 
in direction /i. 

For a generic lattice in d dimensions, the number of 
possible free spectra is thus equal to the number of al- 
lowed topological sectors for the vacuum, i.e., 2**. When 
at least two lattice extensions are equal, the number 
of possible free spectra is reduced. For lattices with 
Lfi = L^ V/i, i^ = 1, . . . , d, the different patterns are la- 
beled by the number of c^ values that arc equal to i, 
thus yielding d + 1 inequivalent possibilities. 

Configurations close to the free limit are expected to 
approach (modulo gauge transformations) one of the pos- 
sible free vacua. This is confirmed by our lattice simula- 
tions, where we find the distribution of P^ to be peaked 
at ±1. The corresponding free vacua can then easily be 
identified by the sign of the averaged traced Polyakov 
loops Pfj, in the various directions. Accordingly, Eq. (|3.2[) 
provides a first approximation to the observed spectrum 
of the staggered Dirac operator. 

This prediction can be refined as follows. For a given 
configuration, let us introduce a configuration built from 
uniform links C/^ (x) = f/^ in each direction, taking values 
in an Abelian subgroup of SU(2) (for instance the diago- 
nal one), and yielding the same averaged traced Polyakov 
loops as the original configuration. Since for uniform and 
commuting links all plaquettes are equal to unity these 
configurations may also be called vacuum configurations. 

For these vacuum configurations, the gauge transfor- 
mation 

Ut,{x) 1-^ g{x) Uf,{x) g\x + fi) 

.1 



operator in such a vacuum configuration is given by 



with 9[x)^\{{U^r 



(3.3) 



can be used in order to trivialize all Polyakov loops to 
P^ = 1 at the expense of introducing periodicity only up 
to g{x) (with Xii = L^). The latter equals the original 
Polyakov loop P^ and (in the diagonal subgroup) com- 
prises two opposite phases, which behave like the con- 
stants c^ (this is a generalization of our previous argu- 
ment that P^ = — 1 can be absorbed by switching c^ 
between and ^). The spectrum of the staggered Dirac 



±i 



\ 



^sir 



27r 



arccos P^ 

"hi 



(3.4) 



Now Cp is again fixed to (i) for (anti)periodic bound- 
ary conditions in the direction /i. Our expectation is 
that, close to the free limit, Eq. (|3.4p provides a better 
approximation than Eq. (j3.2p to the observed spectrum 

of Dks- 

The eigenvalues in Eq. (|3.4p have a multiplicity of 2'*, 
which we derive in the following paragraph. For dimen- 
sion d = 4, this implies an eightfold degeneracy in addi- 
tion to Kramers' degeneracy. A small perturbation of the 
vacuum generically lifts this eightfold degeneracy but not 
Kramers' degeneracy, which remains exact. This mech- 
anism gives rise to what we will call clusters of eight 
eigenvalues below. 

In order to explain how the multiplicity of 2'^ comes 
about, we have to digress to a sketch of the derivation 
of Eq. (|3.4p . The spectrum is most easily obtained by 
looking at the square -D^g of the staggered Dirac operator 
and applying it to plane waves. After some algebra, one 
finds the eigenvalues of D^g to be given by 



A. 



^sin^ 



27r 



arccos Pp 

"hi 



(3.5) 



with < /cp < L^ 



1. Since A„ is invariant under 



kfi ^^ k^ + L^/2, we can restrict the wave numbers to 
< A:^ < Lp/2 — 1 and assign a multiplicity of 2*^ to each 
eigenvalue. An additional multiplicity factor of 2 arises 
from the color degeneracy, yielding an overall multiplic- 
ity of 2^*+^. Because of the symmetry of the spectrum of 
Dks about A = (which arose due to {Dks, 5} = 0, see 
Sec. Ill A 21 above), the eigenvalues of Dks are given by 
the positive and negative square roots of A„ with each 
eigenvalue A„ having half the multiplicity of the corre- 
sponding A„, i.e., 2''. 



B. Adjoint representation 

The situation is similar for the adjoint representation 
of SU(2). Since the latter is insensitive to the group cen- 
ter, all trivial vacua are equivalent to the configuration 
with all links equal to unity. 

The construction above can be repeated by considering 
a vacuum configuration built from link matrices of the 
form 



f/u 



COSOfn 







(3.6) 







where a^ is related to the traced Polyakov loop P^ by 
arccos((Pp - l)/2). The analog of Eq. ^^^ now 



Lfj.a^ 



reads 



A± = ±i 



\ 



T. 



sill 






(3.7) 



where n G { — 1,0,1} and < fc^j < L^/2. In con- 
trast to the situation for the fundamental representation, 
Eq. p.7p predicts that one-third of the eigenvalues re- 
mains unchanged, i.e., they are identical to the eigenval- 
ues in the trivial vacuum (the configuration with all links 
equal to unity). 

An analysis analogous to that for the fundamental rep- 
resentation shows that now the multiplicity is 2''"^ (as 
opposed to 2'^ above). However, here wc have no Kramers 
degeneracy, and thus a small perturbation of the vacuum 
again gives rise to clusters of eight eigenvalues in dimen- 
sion d = 4. 



IV. NUMERICAL RESULTS 
A. Simulation details 

Our numerical results are obtained from sets of 
quenched SU(2) configurations generated using the Wil- 
son gauge action. The simulation algorithm is based on 
a combination of Metropolis and over-relaxation; center 
rotations to explore different topological sectors are im- 
plemented as well. 

We obtain the full spectrum of the staggered Dirac op- 
erator from ensembles of configurations on hyper-cubic, 
isotropic lattices with volumes ranging from V ~ A'^ 
to 16'^. The spectrum is evaluated using the CuUum- 
Willoughby version of Lanczos' algorithm [33]; periodic 
(antiperiodic) boundary conditions are assumed in the 
spatial (temporal) directions. 

For each lattice volume V and for each value of the 
Wilson gauge action parameter /?, our analysis of the full 
staggered spectrum is based on a number of thermalized 
and uncorrelated configurations between a few tens and 
a few thousands. For each configuration the number of 
distinct eigenvalues with positive imaginary part is V/2 
for the fundamental representation and 'SV/2 for the ad- 
joint representation. Hence, our data for each pair (V, /3) 
typically contains on the order of one million distinct 
eigenvalues. 

Furthermore, we also generate gauge configurations on 
much larger lattices (up to 34 x 38 x 46 x 58) in order to 
investigate the distribution of eigenvalues from Eqs. p.4p 
and (|3.7[) . On these lattice we do not diagonalize the 
Dirac operator but calculate only the averaged Polyakov 
loops. 



B. Separation of spectral scales 

When f3 is increased to large values, on a lattice with 
a fixed number of sites in each of the four directions, 



the spectrum of Dks shows structure on three different 
scales, see Figs. [T] and [5] for examples with fundamental 
and adjoint fermions, respectively. 



1 . Plateaux 

The coarse structure is given by the free limit, i.e., 
the eigenvalues approach the values of Eq. (|3.2p . with 
c^ £ {O, 5} chosen to reflect the sign of the Polyakov 
loops as described below Eq. (|3.2[) . The collection of 
all eigenvalues in the vicinity of the values predicted by 
Eq. p.2p we call a plateau, indicated by dashed blue lines 
in Figs. [Hand H 

As discussed in Sec. IHI Al for the fundamental repre- 
sentation there are five different plateau structures. We 
find all these classes in our simulated configurations and 
show representatives for each class in Fig. [TJ (The spec- 
tra shown are the first representatives of each class which 
we come across in the Monte Carlo history.) 

For adjoint fermions the plateau structure is uniquely 
determined by the lattice sizes (and the boundary con- 
ditions) since the adjoint representation is center-blind, 
as discussed in Sec. IHIBI Thus, configurations with dif- 
ferent signs of P^ lead to the same plateau structure, see 
Fig. [3 



2. Clusters 

A closer look at the staggered spectra reveals that the 
distribution of eigenvalues inside a given plateau shows 
additional structure. The eigenvalues are grouped in 
clusters of eight. When the free limit is approached, 
the typical separation between nearest clusters within the 
same plateau is smaller than the separation between dif- 
ferent plateaux, but larger than the separation of eigen- 
values within each cluster. 

The position of the clusters can be predicted by 
Eqs. (|3.4p and (|3.7p for the fundamental and adjoint rep- 
resentation, respectively. In Figs. [T] and [2| these predic- 
tions are indicated by solid red lines. Thus, the clus- 
ter structure is determined by the traces of the averaged 
Polyakov loops (and the lattice size and boundary con- 
ditions). For the adjoint representation, we see that one- 
third of the eigenvalues forming a plateau does not split 
into clusters, but stays close to the plateau levels. This 
is also in agreement with our discussion of Eq. (|3.7p . 

Recall that, approaching the free limit, the distribu- 
tion of the traced Polyakov loops becomes peaked at ±1, 
corresponding to the center elements of SU(2). This a 
posteriori justifies the approximation of the staggered 
spectrum by the plateaux. 

The observation that each cluster contains eight eigen- 
values reflects the fact that the eigenvalues of Eqs. (|3.4p 
and (|3.7p are degenerate with multiplicity 8. These eigen- 
values where calculated for vacuum configurations with 
uniform and commuting links t/p(a;) in each direction. 
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FIG. 1: (Color online) Eigenvalues of Dks close to the free limit for the fundamental representation of SU(2) (L** = 4^*, /3 — 
1000). We display representatives for the five different plateau patterns (indicated by dashed blue lines) described below 
Eq. p.2|) . Within the plateaux the eigenvalues (black dots) arrange themselves in clusters of eight, whose locations are predicted 
by Eq. p.4|) (solid red lines). The lower right panel shows a spectrum from a 6'^ lattice, with a richer cluster structure. 
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FIG. 2: (Color online) Eigenvalues of Dks close to the free limit for the adjoint representation of SU(2) {L = 4 , /3 = 1500). 
The unique plateau structure is marked by dashed blue lines. Within the plateaux the eigenvalues (black dots) arrange 
themselves in clusters of eight, whose locations are predicted by Eq. (|3.7|l (solid red lines). 



Since the simulated configurations close to the free limit, 
which we approximate by these vacuum configurations, 
do not have exactly uniform and commuting links, the 
eightfold degeneracy is lifted. 

In the lower right panel of Fig. [T] we also show an ex- 
ample obtained from simulations on a 6^ lattice. We see 
that for larger lattices the expected patterns get more 
and more complicated, but the agreement persists with 
just the four parameters P^ of Eq. p.4p determining the 



complete cluster structure. 

C. Spectral statistics 

1. Unfolding 

Before one can discuss spectral statistics and compare 
to. e.g.. predictions from RMT, the spectra have to be 
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FIG. 3: (Color online) Level spacing densities obtained from spectra of the staggered Dirac operator in the fundamental 
representation of SU(2) for small j3. The numerical data are consistent with the Wigner surmise for the chSE - as expected 
according to the symmetry properties of the staggered Dirac operator, see Sec. Ill Al 



unfolded, i.e., the (imaginary parts of the) eigenvalues 
have to be rescaled such that the mean separation be- 
tween adjacent levels is unity. This can be achieved by 
defining a new spectrum {a;„} with a;„ := N{Xn), where 
A^(A) denotes the mean integrated spectral density, i.e., 
A^(A) is a smooth function satisfying A^(A) w #{A„ < A} 
in some approximate or asymptotic sense. 

In low-dimensional quantum chaos, one usually has an 
analytical formula for N{X), provided by Weyl's law, see, 
e.g., [23]. For the QCD Dirac operator first steps towards 
an equivalent asymptotics have been discussed in pJ8| . 
Until now, however, this approach does not provide a 
formula which can be used for unfolding spectra in lattice 
QCD. 

If one has no analytical prediction for the mean spec- 
tral density, it has to be extracted from the data them- 
selves. The latter can be done by averaging over sev- 
eral spectra, which is known as ensemble unfolding, thus 
yielding one mean density for a whole ensemble of spec- 
tra. In contrast, fitting an ansatz to the spectral density 
of an individual spectrum or extracting its mean density 
by a moving average in A is known as configuration un- 
folding. The resulting mean densities will in general differ 
from each other. This ambiguity can make it difficult to 
extract reliable information on long-range spectral corre- 
lations, see the detailed discussion of different unfolding 
methods for lattice Dirac spectra in Ref. [3^ . The statis- 
tical function we are interested in is always the density 
P{s) of nearest-neighbor spacings. Since P(s) measures 
short-range spectral correlations, changing the unfolding 
method will not impair our results, as long as we make 
sure that each unfolding procedure is stable and consis- 
tent. 

As we want to discuss spectral statistics for very dif- 
ferent values of the gauge action parameter f3 and on dif- 
ferent scales, we are forced to employ different unfolding 
methods, each one of them tailored to fit the require- 
ments of the particular situation. 

For large /3 the spectra show different scales, see 
Sec. IIVBI Therefore, we unfold separately on each scale. 



The different plateau and cluster structures forbid meth- 
ods of ensemble unfolding, and thus we employ configu- 
ration unfolding. The methods we use are all variants of 
what is called local unfolding in Ref. [39|. 

When studying spacings within clusters we divide spac- 
ings of adjacent levels by the mean level spacing within 
the cluster. The latter we calculate as the difference be- 
tween the largest and the smallest level divided by seven. 
The unfolded spacings then have unit mean as required. 

For spectral statistics between clusters within a pla- 
teau, we proceed analogously. We divide spacings be- 
tween adjacent clusters by the difference between the po- 
sition of the largest and the smallest cluster divided by 
the number of clusters within the plateau minus one. 

Finally, when we want to study spacings between pla- 
teaux, we adapt the previous methods as follows. Wc di- 
vide the difference between the position of plateau j + 1 
and plateau j by the difference between the positions of 
plateaux j + 5 and j — 4 divided by the number of plateau 
spacings in this range, i.e., by nine. 

For small /3 the spectra do not show different scales. 
Thus, complete spectra could be unfolded in one go, and 
both configuration and ensemble unfolding are admissi- 
ble. We experimented with the ensemble unfolding de- 
scribed in Ref. [I^l, which works reliably for small /3. 
However, when we increase (3 the method has to break 
down, due to the different plateau structures which begin 
to emerge. In order to have an unfolding method which 
does not exclude an intermediate /3-range but instead al- 
lows us to go smoothly from small to large /3, we decided 
to unfold spectra for small /3 in the same way as we un- 
fold spacings within clusters for large /3. Effectively, this 
means that we neglect every eighth spacing and unfold 
the remaining spacings on the scale of the neighboring 
eight eigenvalues. 



2. Level spacings for small f3 

For small values of /?, i.e., far away from the free and 
continuum limits, the level spacing densities have to be 
compared to the predictions reported in Sees. Ill A~2l and 

lnB2l 

Figure [3] shows the level spacing densities obtained 
from spectra of the staggered Dirac operator in the funda- 
mental representation of SU(2). The numerical results on 
all our lattices are consistent with the chSE, in contrast 
to the chOE, which would be expected in the continuum. 

Likewise, in Fig. [4] we show level spacing densities ob- 
tained from staggered Dirac spectra in the adjoint rep- 
resentation of SU(2), which now are consistent with the 
chOE as expected. Again this is in contrast to the contin- 
uum situation in which one should have chSE statistics. 

This behavior of the staggered Dirac operator in SU(2) 
gauge fields, having spectra in different universality 
classes than the corresponding continuum operator, has 
been observed earlier [l2|, |14|, ilT,] . 



3. Level spacings for large j3: Approaching the free limit 

For large /3, plateaux and clusters emerge. Therefore, 
we discuss spectral statistics separately on the different 
scales. 

(a) Spacings within clusters. Generically, the eightfold 
degeneracy of the levels predicted by Eqs. p.4p and (|3.7p 
is lifted by the nonuniformity of the gauge field configu- 
ration. For this reason we observe small clusters of eight 
eigenvalues. 

Close to the free limit we can think of treating the 
nonuniformity of the gauge field configuration as a small 
perturbation of the corresponding vacuum configuration, 
i.e., of the configuration with uniform and commuting 
links leading to the same Polyakov loops. A cluster then 
arises by diagonalizing an 8 x 8 matrix, the perturbation 
restricted to the subspace corresponding to a degener- 
ate eigenvalue of Eq. p.4p or Eq. p.7p . This matrix 
inherits symmetries and effective randomness of the per- 
turbation, i.e., of the gauge field part of the staggered 
Dirac operator. Therefore, we expect the distribution of 
level spacings to follow the same chRMT prediction as 
for small /?. 

Our data for fermions in the fundamental representa- 
tion of SU(2) confirm this expectation. As Fig. [5] shows, 
the level spacing density within each cluster is consistent 
with the prediction from the chSE. Note that this is true 
over a very large range of /3 values. 

Also in the adjoint representation the spacings within 
clusters follow the same pattern as the spacings for small 
/?, in this case leading to a chOE distribution, see Fig. [B] 
Note that for the adjoint representation, the clusters at 
the ends of a plateau show some nongeneric features 
which would require a more sophisticated unfolding pro- 
cedure. We avoid this complication by restricting the 



analysis to cluster spacings in the central part of each 
plateau. 

These results confirm once more that the spectral prop- 
erties of the staggered Dirac operator are different from 
those of the continuum operator. Moreover, in the sense 
described above, our analysis demonstrates that this dis- 
crepancy persists when approaching the free limit. 

(b) Spacings between clusters. For the fundamental 
representation the cluster positions are, to a good ap- 
proximation, described by Eq. (|3.4p . Both Eq. (|3.2p and 
Eq. (|3.4p describe spectra of integrable systems, and thus 
one generically expects uncorrelated levels, like for a Pois- 
son process [2J, [2^ ■ 

However, on commensurate lattices, i.e., on lattices 
with rationally dependent i^/2, Eq. (|3.2[) predicts many 
accidental degeneracies which would lead to nongeneric 
spectral statistics. Generic behavior can in principle be 
restored in two different ways: On the one hand by chang- 
ing to lattices with rationally independent extensions 
i^/2, and on the other hand by introducing additional 
phase shifts, like the Polyakov loops do in Eq. (j3.4p . How- 
ever, as we pointed out earlier, close to the free limit the 
distribution of averaged traced Polyakov loops is peaked 
at ±1. In Fig. 3 of Ref. [2y] it has been observed that 
these almost equal phase shifts in each direction are not 
able to restore generic behavior. 

Therefore, we choose to study a large incommensurate 
lattice with extensions Li x L2 x ^3 x ^4 ~ 34 x 38 x 
46 X 58. We refrain from diagonalizing the Dirac oper- 
ator on this lattice but instead calculate the averaged 
traced Polyakov loops and determine the approximate 
cluster spectrum from Eq. (|3.4p . Then the density of 
spacings between different clusters within the same pla- 
teau agrees well with the prediction from a Poisson pro- 
cess, Ppoisson(s) = e"'*, see Fig. [7] (left). The same holds 
true for fermions in the adjoint representation with the 
approximate cluster spectrum determined from Eq. (|3.7p , 
see Fig. [3 (middle). 

Note that we analyze the spacing distribution for the 
spectra p.4p and (|3.7p but not the spacings between the 
actual cluster positions which one would obtain by av- 
eraging over all eigenvalues belonging to a given cluster. 
As can be seen in Figs. [1] and [2] these actual positions 
differ slightly from the predicted values. Whether or not 
these differences could lead to a deviation from Poisson 
behavior is an open question. 

(c) Spacings between plateaux. The positions of the pla- 
teaux can be approximately described by the free spec- 
trum, Eq. (|3.2p . As in the case of spacings between clus- 
ters, which we discussed above, we thus once more study 
the spectrum of an integrable system. Again accidental 
degeneracies lead to nongeneric statistics on commensu- 
rate lattices. Generic behavior can be restored by switch- 
ing to an incommensurate lattice, which we demonstrate 
in Fig. [7] (right). Level spacing distributions like for a 
Poisson process were predicted and observed earlier for 
incommensurate lattices [40, |41| . 

We add the same note of caution as for the spacings 



l" = 8"*, P = 1.5, 100 configurations 



L = 8 , P = 5, 100 configurations 




FIG. 4: (Color online) Level spacing densities obtained from spectra of the staggered Dirac operator in the adjoint representation 
of SU(2) for small /3. The numerical data are consistent with the Wigner surmise for the chOE - as expected according to the 
symmetry properties of the staggered Dirac operator, see Sec. IIIBI 
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FIG. 5: (Color online) Level spacing densities for eigenvalues within clusters, obtained from spectra of the staggered Dirac 
operator in the fundamental representation of SU(2) for increasing p. Agreement with the chSE persists for large /3, i.e., close 
to the free limit. 



between clusters. Our analysis concerns the level spac- 
ings of the spectrum (|3.2p but not the spacings between 
the actual plateau positions, obtained by averaging over 
all eigenvalues within a given plateau. 

Let us also point out that, in spite of the strong simi- 
larities between the data describing the cluster spacings 
and the plateau spacings in Fig. [71 they describe correla- 
tions on spectral scales that typically differ by an order 
of magnitude. 



V. CONCLUSIONS AND OUTLOOK 

We have studied the spectral properties of the stag- 
gered Dirac operator Dks when approaching the free 
limit for gauge group SU(2), both in the fundamental 
and in the adjoint representation. With SU(2) gauge 
fields the staggered Dirac operator on the lattice belongs 
to a different symmetry class than the corresponding con- 
tinuum operator. 

Our numerical analysis revealed that, when the free 
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L* = g'', p = 50, 100 configurations 



L* = 8", p = 150, 100 configurations 
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FIG. 6: (Color online) Level spacing densities for eigenvalues within clusters, obtained from spectra of the staggered Dirac 
operator in the adjoint representation of SU(2) for increasing /3. In the central part of each plateau (cf. text) agreement with 
the chOE persists for large /3, i.e., close to the free limit. 
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FIG. 7; (Color online) Left; The spacing density of the approximate cluster spectrum of the fundamental representation 
predicted by Eq. (|3.4p for a single configuration on a 34 x 38 x 46 x 58 lattice at /3 = 10000 agrees with the spacing density 
e'" for a Poisson process. Middle: Same as left panel, but now for the adjoint representation and Eq. (|3.7|l . Right: The 
spacing density of the plateau positions, i.e., the eigenvalues of the free staggered Dirac operator, as predicted by Eq. (|3.2p on 
a 34 X 38 X 46 X 58 lattice also agrees with the spacing density e~" for a Poisson process. 



limit is approached, the spectrum of the staggered Dirac 
operator shows structure on three well-separated scales. 
The behavior on the two coarser scales is characterized 
by the formation of plateaux and clusters. 

We have shown that for a given gauge field configura- 
tion the positions of plateaux and clusters can be pre- 
dicted analytically. To this end we have constructed a 
vacuum configuration with uniform and commuting links 
in each direction, chosen such that they reproduce the 
averaged traced Polyakov loops of the original configura- 
tion. The Dirac spectrum in this vacuum configuration 
yields a good approximation to the plateau and cluster 
positions. In turn, plateaux and clusters are determined 
by lattice geometry, boundary conditions, and Polyakov 
loops alone. 

Our model also predicts a systematic degeneracy of the 
staggered spectra in nontrivial vacuum configurations: 
All eigenvalues have a multiplicity of eight (in addition to 
Kramers' degeneracy in the fundamental representation). 
For typical gauge field configurations this degeneracy is 
lifted leading to the formation of clusters of eight eigen- 



values, which we observe numerically. 

We have analyzed spectral correlations on all three 
scales in terms of the distribution of spacings between 
adjacent eigenvalues, clusters, and plateaux. Spacings 
between approximate plateau and cluster positions are 
uncorrelated as for a Poisson process, whereas level spac- 
ings on the finest scale, i.e., within clusters, follow the 
chRMT predictions. For the latter the symmetry class is 
always that of the staggered operator and never that of 
the continuum operator, even for very large /?, i.e., close 
to the free limit. 

Finally, we briefly comment on the possibility to re- 
cover the symmetry class of the continuum operator in 
spectra of the staggered Dirac operator in SU(2) gauge 
fields. We restrict ourselves to the fundamental represen- 
tation. For every finite lattice spacing, the antiunitary 
symmetries of the i'KS operator are different from those 
of the continuum Dirac operator. When approaching the 
continuum limit, i.e., when increasing (3 at fixed physi- 
cal volume, the eigenvalues should form near-degenerate 
quartets (more precisely, near-degenerate pairs of exactly 
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Kramcrs-dcgcncratc pairs), signaling the suppression of 
imphysical taste-changing interactions of the staggered 
operator. This has been seen for the first time in Ref. [23| 
using highly improved staggered fermions. In |23| it was 
also shown that the distribution of the smallest eigen- 
value makes a transition from chSE to chOE, i.e., from 
the symmetry class of the lattice operator to that of the 
continuum operator. Related phenomena have been ob- 
served for improved staggered fermions in the fundamen- 
tal representation of SU(3) [H, El, Ell . We expect that 
the distribution of spacings between eigenvalues will also 
show a similar transition to the continuum behavior. Our 
present findings do not contradict such an expectation as 
our study corresponds to a different physical setting: In- 



stead of the continuum limit we have investigated the 
behavior in the free limit, where RMT applies only at 
the finest scale of eigenvalue fluctuations. 
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